Charge Fluctuations and the Valence Transition in Yb under Pressure 
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We present a dynamical mean field theory study of the valence transition {f^^ /^^) in elemental, 
metallic Yb under pressure. Our calculations reproduce the observed valence transition as reflected 
in the volume dependence of the 4/ occupation. The transition is accelerated by heating, and 
suggests quasiparticle or Kondo-like structure in the spectra of the trivalent end state, consistent 
with the early lanthanides. Results for the local charge fluctuations and susceptibility, however, 
show novel signatures uniquely associated with the valence transition itself, indicating that Yb is 
a fluctuating valence material in contrast to the intermediate valence behavior seen in the early 
trivalent lanthanides Ce, Pr, and Nd. 
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PACS numbers: 71.20.Eh, 74.40.+k, 75.20.Hr 



The valence state of rare earth atoms in lanthanide 
compounds has a crucial effect on their physical proper- 
ties. Determination of the lanthanide valence from first 
principles and description of the 4/ electrons has been a 
long standing challenge due to the duality between their 
atomic character, stemming from the on-site electron- 
electron interaction, and the itinerant character, due 
to the lattice periodicity. Theories based on the two 
species picture, which treat part of the 4/ electrons as 
atomic and the rest as itinerant, succeeded in reproduc- 
ing the trends across the lanthanide series for compounds 
with integer valence [J] . Nevertheless, besides being con- 
ceptually unsatisfactory, the two species picture cannot 
describe transitions between different valence states as 
well as such outstanding behavior as heavy fermions. 
Only recently has the combination (LDA+DMFT) [3,[ll 
of the local-density approximation (LDA) and dynami- 
cal mean-field theory (DMFT) been applied to the rare 
earth materials [1, d, H, 0, d, 0| , providing a material spe- 
cific electronic structure including the local many-body 
dynamics. 

In this Letter we employ the LDA-I-DMFT method to 
study the valence transition and charge fluctuations in 
elemental Yb metal. It is well known that Yb and Eu 
behave differently from the other lanthanides in their ele- 
mental form. If we define the valence state as the number 
of electrons participating in bonding, the majority of the 
lanthanide series is trivalent, however for Yb and Eu the 
3+ and 2+ valence states are close to degenerate with 
2+ state being more stable at ambient conditions [ij. 
This results in a number of anomalous properties, such 
as a larger molar volume (as compared to the trend the 
rest of the lanthanide series follows) , a lower bulk mod- 
ulus [l^, and the thermal expansion coefficient of Yb 
being three times larger than for most other lanthanides 
The gradual transition from the divalent to the 
trivalent state occurs in Yb over the range 0-34 GPa, 



where the full 4/ shell is opened as an / electron is pro- 
moted to the valence band. The doped 4/ holes can move 
through the crystal by thermally activated hopping with 
spd bands acting as particle and energy reservoirs, fluc- 
tuating valence behavior, or they can move coherently 
between atomic sites forming a narrow band well known 
for heavy fermions, intermediate valence behavior. We 
will discuss how these concepts apply to Yb. 

In the present study we start with a self-consistent 
LDA calculation, transform its one-electron Hamiltonian 
into an orthogonal Wannicr function basis, H^^^, and 
then calculate the 4/ interaction parameter U as well 
as the the double counting correction to H^^^. Since 
the valence transition of Yb is of principal interest here, 
we have to retain the 6s, 6p, and 5d valence orbitals to 
allow changes of the /-shell occupation. We assume the 
SU(N) symmetric form of the local interaction 
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(the interaction parameter U is a. function of the spe- 
cific volume is shown in the inset of Fig. [1]), and express 
the local Green's function in the relativistic j = 5/2, 
7/2 basis to make the off-diagonal elements small so 
that they can be neglected. To solve the auxiliary im- 
purity problem we employ two Quantum Monte Carlo 
(QMC) solvers, one using the Hirsch-Fye (HFQMC) ^ 
algorithm and the other using the hybridization expan- 
sion continuous time QMC algorithm (CTQMC) [13[ as 
well as an approximate, but computationally efficient 
Hubbard-I (HI) solver. Extrapolations of the HFQMC 
results to an infinite number of imaginary time slices 
L were found to agree within statistical uncertainties 
with the CTQMC. However, we point out that espe- 
cially at the larger L values, the HFQMC calculations 
had significant difficulties with ergodicity in the midst 
of the Yb valence transition where large fluctuations in 
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FIG. 1: (Color online) The / shell occupation vs pressure 
at various temperatures compared to the experimental data 
(closed symbols). Open symbols are QMC data; the 630 K 
and 1160 K curves were obtained with CTQMC, and the 
1580 K curve with HFQMC at L = 80. The consistency 
of HFQMC with the CTQMC data was checked by L ^ oo 
extrapolation of the L =80, 112 and 160 data at 630 K. X-ray 
absorption edge spectroscopy (XAES) data was taken from 
Ref. [l^, and resonant inelastic x-ray scattering (RIXS) data 
taken from Ref. [l^. Inset shows the interaction U vs. specific 
volume. 



71/ were encountered in the Ising space sampling. The 
CTQMC method did not encounter such problems, nor 
did HFQMC without any valence transition (e.g., Ce, 
Pr, and Nd 0,[ii). 

In Fig. [1] we show the evolution of the /-orbital occu- 
pation n/ as a function of the specific volume V at var- 
ious temperatures. Comparison to the experimental x- 
ray absorption edge spectrosocpy (XAES) and resonant 
inelastic x-ray spectroscopy (RIXS) data shows some un- 
derestimation of Uf at lower volumes, possibly originat- 
ing from the lack of charge self-consistency, which would 
take into account the correlation driven redistribution 
of charge between the / and spd orbitals. Overall the 
DMFT(QMC) results agree reasonably well with the ex- 
perimental data over the entire pressure range of interest. 
Note that the rising 6s and 6p bands lead to an increase 
in n/ with compression for the trivalcnt rare earths 0, 
so that it is hkcly that the plateau in the DMFT(QMC) 
value reached at the smallest volumes in Fig. [T] is indeed 
near the end of the transition from divalent to trivalcnt 
in spite of the fact that is still larger than 13.0. The 
ability of DMFT(HI) to produce the valence transition 
seen in experiment suggests that the physics is essen- 
tially related to charge transfer between the / and spd 
orbitals. At smaller volumes where delocalization of /- 
electrons becomes more important HI overestimates the 
number of holes in the 4/ shells. 

The calculated Uf values show a sizable temperature 
dependence both with HI and QMC, particularly at low 
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FIG. 2: (Color online) Plot of Yb spectral function A{uj) 
calculated in (a) HI at T = 630 K and (b) CTQMC at T = 
1160 K for selected volumes (in A'' /atom). Dashed curves 
in (b) are for i = f ; solid curves, j = |- Spin-orbit split 
peaks corresponding to the f^'^ — > f^'^ excitation are near the 
chemical potential, while those corresponding to the /^^ — > 
/^^ excitation are around —U{V). An expanded view near 
the Fermi-level for V = 20.3 A'^/atom is shown in (c). The 
j = I curve shows two shoulders around the large central 
peak. 



pressure as is evident in Fig. [TJ The general trend is that 
at higher temperatures, the f^^ state is favored over the 
f^^ state. This trend can be followed down to the low- 
est studied temperature of 16 K, which is only accessible 
with HI. We expect this trend would hold at lower tem- 
peratures for QMC calculations as well. Temperature de- 
pendent measurements of the valence of Yb in YbInCu4 
show the same trend (4/ occupation decreasing with in- 
creasing temperature) [l7[ , with different studies finding 
a minimum of about 13.1 [isl . [l9| . The same behavior 
was found for YbAgCu4 and YbAls . The authors of 
Ref. did not find any significant difference in elemen- 
tal Yb between the temperatures of 250 K and 25 K; at 
250 K Yb should be already strongly divalent and low- 
ering the temperature cannot decrease the valence any 
further. There is, however, a significant broadening of 
the spectrum at higher temperature, so we expect that 
the trend seen in our calculations would be borne out if 
measurements were carried out at higher temperatures. 
At high pressure the temperature dependence of Uf be- 
comes weak. The chemical environment of Yb can affect 
the temperature sensitivity of Uf as well; in fact, the 
measured valence of Yb in YbCaGe has been recently 
found to be temperature independent plj , and Yb nearly 
divalent. This is somewhat a contrast to our result here, 
where the temperature sensitivity is only near the di- 
valent state. The details of the electronic structure of 
YbCaGe that causes the Yb valence to be temperature 
independent are not understood. 

The calculated 4/ spectral densities A{lli) are shown 
in Fig. [21 as obtained by direct evaluation of the Creen's 




FIG. 3: (Color online) (a) Equal time charge fluctuations in 
Yb (CTQMC), Ce, Pr and Nd (HFQMC with L ^ oo extrap- 
olation) at 630 K. If not shown error bars are smaller than 
the symbol size, (b) Charge fluctuations for Yb (diamonds, 
the same as in left panel) compared to their imaginary time 
average, Tx^-^\ for Yb (closed circles) and Ce (open circles). 



function near the real frequency axis for HI, and applying 
the maximum entropy method [S^l to the CTQMC data. 
At low pressure, there is a clear spin-orbit-spht peak near 
the chemical potential (w = 0), which corresponds to the 
/^^ f^^ excitation. As the pressure is increased, the 
weight of this double-peak decreases, and the /^"^ /^^ 
excitation appears around 8 eV (w U) below the f^'^ 
f^"^ double-peak. While the analytic continuation smears 
the high energy features of the CTQMC spectra. HI is 
known to incorrectly reduce the width of the Hubbard 
bands. 

When the valence transition approaches completion 
as in the HI result (bottom of Fig. [2^), the j = 7/2 
/^^ f^^ peak overlapping the Fermi level shifts in- 
creasingly above this energy to become the unoccupied 
upper Hubbard band for the 4/ states, corresponding to 
the new 4/ hole. A close examination of the CTQMC 
counterpart (Fig.[2t), on the other hand, shows the cor- 
responding j = 7/2 structure to contain a relatively 
large peak at the Fermi level with shoulders to either 
side. Based on the temperature scaling of the Fermi- 
level peak, we believe it to be a quasiparticle or Kondo- 
like contribution consistent with increasing delocaliza- 
tion of the 4/ hole as volume is reduced, and that 
the lower shoulders are likely of more Hubbard char- 
acter. Moreover, any transfer of spectral weight into 
such a quasiparticle contribution at the expense of the 
/^^ — > f^^ Hubbard band would also contribute to a 
larger DMFT(QMC) n/ than in DMFT(HI) as is ob- 
served. Finally, both the simultaneous existence of two 
lower Hubbard bands, and the gradual transfer of spec- 
tral weight evident in Fig. [2] are fundamental signatures 
of electron correlation which a single-particle approach 
such as LDA or LDA-I-U cannot reproduce. 



Thus far we have shown a gradual emptying of th§ 
4/ shell with applied pressure and the correspond- 
ing changes in the single-particle spectra interpreted in 
terms of the relative abundance of f^^ and /^^ config- 
urations in the system. We now present an analysis of 
the associated local charge fluctuations in Yb. In Fig. [3^ 
we compare the local charge fluctuations in Yb, defined 
through the mean square deviation of the /-occupation 
{Snj); Sfif = fif — Uf , with those in Ce, Pr and Nd [3]. 
While Ce, Pr and Nd show a monotonic increase of (Jn^) 
with pressure, the charge fluctuations in Yb exhibit a 
local maximum around 27 A'^/atom where ~ 13.5. 
Evaluating the thermal expectation in a basis of n/ 
eigenstates yields {5nj) = u;i3(l — W13), if we assume 
only f^^ and /^^ configurations have significant weight 
(ifi3 + wi4 = 1), and thus a peak in the midst of the 
valence transition 71113 = — *■ 1. On the other hand, 
if a single /" configuration predominates (w„ '^1) but 
then loses weight to /"^^ and /"+^ states due to growing 
hybridization under compression, one will see a mono- 
tonic increase in (Snj:). It would appear likely that the 
behavior seen in Fig. [3^ for Yb is reflects a combination 
of both effects, while that of Ce, Pr, and Nd is primarily 
the latter delocalization behavior. 

To gain more insight we further evaluate the local 
charge susceptibility x*"^"* (Fig. [Sb), obtained from the 
imaginary-time correlation function 

X^'^ = r dr x{t), x{r) = {Shf{T)Shf{0)). (2) 
Jo 

We find that the initial increase of the charge fiuctua- 
tions with pressure is mirrored by an increase of the local 
charge susceptibility with both having their maxima at 
about the same volume, however, the increase of (Snj:) 
at high pressure is not reflected in the charge suscepti- 
bility. The origin of this behavior is revealed in the inset 
of Fig. m where we compare the imaginary time charge 
correlations for 1^=13.6 and 29.0 A^/atom. While the 
magnitude of the charge fluctuations is about the same 
(r = intercepts); the fluctuations at high pressure are 
short lived (rapid decay with increasing t) leading to a 
relatively lower susceptibility. In Fig.|3]we show the cor- 
responding physical susceptibilities on the real frequency 
axis, imaginary parts of which characterizes the density 
of charge excitations. 

Xir) = - / dco x^^H^) (3) 

TT Jo 1 - e 

Comparison of x(t) and x*-^-* (ijj) shows that slowly decay- 
ing fluctuations are, as expected, related to low energy 
charge excitations. 

Low energy local charge excitations in strongly corre- 
lated metals arc rather rare since they arc usually sup- 
pressed by the onsitc Coulomb repulsion, as is the case 
for Ce, Pr, and Nd where Tx^^-* is 20-60 times smaller 
than {Snj:) over the volume range of Fig. [4)d [i3|. Exis- 
tence of such excitations and the corresponding maxima 
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FIG. 4: (Color online) Imaginary part of the local // charge 
susceptibility at 630 K and specific volumes of 13.6 (black 
dashed line) and 29.0 (red solid line) per atom. The inset 
shows the corresponding imaginary time correlation functions 
obtained with CTQMC. 



in x^^-* vs V and ((5nj) vs V curves distinguishes this fluc- 
tuating valence behavior of Yb from the intermediate va- 
lence behavior of the other studied lanthanides. The cal- 
culated local charge susceptibility provides a clear phys- 
ical meaning to the earlier reports of near degeneracy 
of 2+ and 3-|- valence states in Yb There is also a 
direct connection to the experimental observation of the 
unusual pressure-volume curve of Yb, namely the com- 
pressibility of the electronic subsystem is proportional 
to its charge susceptibility. It is plausible to expect that 
the maximum of local x'^"* vs V leads to softening of 
the entire electronic liquid, although the effect will be to 
some extent masked by the presence of spd bands. 

Our local charge fluctuation and susceptibility results 
suggest that pressure affects metallic Yb essentially in 
two ways which can, to a first approximation, be looked 
at separately: hole doping of the filled 4/ shell and grow- 
ing hybridization of the 4/ bands with the spd bands. 
In regard to the first mechanism, it is the crossover of 
atomic-like f^^ and f^^ levels which drives the valence 
transition creating a peak in {Snj^) as noted. Since the 
system is then also sensitive to external perturbations 
breaking this near degeneracy, a large charge suscepti- 
bility also follows. Hole doping into the atomic-like 4/ 
shell can therefore explain both the (local) maxima in 
{Sn^f) and x^^-*. The second mechanism, growing hy- 
bridization with its concomitant 4/ delocalization, leads 
to an increase of charge fluctuation {5nj) at higher pres- 
sures, as well as to the screening of these fluctuations, 
and thus the small charge susceptibility x^^\ in anal- 
ogy to screening of magnetic moment in the Anderson 
impurity model. 

We have examined the pressure induced valence tran- 



sition of elemental Yb using the LDA-I-DMFT approacl^ 
The transition is advanced also by increasing tempera- 
ture at larger volumes, and appears to have reached its 
trivalent limit by the smallest volumes considered, where 
there is evidence of Kondo-like structure in the spectra as 
found with the early lanthanides. Perhaps most interest- 
ingly we find that while the / orbital occupation and the 
single-particle spectra evolve monotonically with com- 
pression, the charge fluctuations and local charge sus- 
ceptibility exhibit a distinct peak. We interpret this fea- 
ture as fluctuating valence behavior in contrast to Ce, Pr 
and Nd which can be classified as intermediate valence 
systems. 
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